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ABSTRACT 


A theoretical model is developed for the wind velocity 
field in the atmospheric near-water layer above ocean waves. 
A turbulent exchange coefficient due to the wave-caused 
Me eeee fluctuations is introduced into the equations to 
account for the corresponding turbulent nature of the motion 
and the resulting energy fluxes. The equations of motion 
for air are derived in a manner to allow for variation of 
the turbulent exchange coefficient with height. <A stream 
function is introduced and a numerical solution obtained by 
utilization of the Richtmyer method. The profiles of wave- 
caused Reynolds stresses and their resulting energy trans- 


fers are discussed and compared with observational data. 
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IL. INTRODUCTION 


Energy and momentum transfer to water waves are pea 
marily dependent upon the conditions of the atmosphere 
immediately adjacent to. the ocean surface. In the case 
of air flow above an open body of water the atmosphere- 
ocean interactions give rise to wave perturbations in the 
upper layer of the water and in the adjacent atmosphere. 
It has been suggested that the wave perturbation charac- 
teristics be used to determine the energy and momentum 
transfers across the air-sea boundary. 

Investigation of the nature of these energy and momen- 
tum transfers have been the purpose of numerous studies 
in recent years. Such studies of wind-wave Coupdiaceweke 
initiated with the quasi-laminar model of Miles (195/94 
who allowed the wave-induced motion to interact with the 
shear flow only at the critical level. Models progressed 
through the studies of Phillips (1957) and Benjamin (1958) 
fo those more recently of Yefimov (1970) and Reynolds 
(1972), wherein there is also interaction with the turbu- 
lent Reynolds stresses. 

In the present study, wave related momentum transfers 
were investigated through the use of a numerical model. 
This study is similar to that of Yefimov (1970) jeep dita ls 
more general and allows for dependence of the mean wind 


on time and variation of the turbulent coefficient WLth 


10 





height. This model utilizes an initial value numerical 
procedure similar to that of Newman (1969). 

In addition, specific data derived from the statisti- 
cal study of Davidson and Frank (1973) will be utilized 
to define solution parameters to generate the numerically 


derived stress patterns. 


UAL 





ETI, FORMULATION OF THE MODEL 


The basis of the problem is the equation of motion. 
Since the stress is expected to be independent of y, two- 
dimensional motion is assumed with the resulting equation 


me che f£orm: 





3° + We vWe - 2 w+ FF, (aa) 
where 

FE S— (ah + a (u'w'), and 

Fo = SS (-utw) + 2 (ut?) 


Here x is the horizontal axis; zis the vertical axis; 0 
moethe density of air; P ws the non-hydrostatic pressure; 
u’' and w' are the turbulent velocity fluctuations. The 
stationary part of each variable (time average) will be 
denoted with a straight line above it and the nonstation- 
ary part (space average) by a waved line. The stationary 
velocities do not include the random components which 
would disappear with a suitable time average. The veloc- 


ity (WV) will be defined as: 


V= (utU)h®+twhk (2.2) 
, / 
where U = U(z) and the factor = will be neglected in the 
development of the equations. For simplicity the problem 


will be restricted to those cases of a single wave number. 
This neglects the effect of wave self-interaction which 


would affect other waves, 


lez 








The turbulent velocity fluctuations, u' and w' are 
assumed to have a dual nature. They are generated by the 
gradient of the mean velocity in the near water layer and 
also by the wave-caused perturbations of the velocity. 
The tangential Reynolds stress (u'w') will be expressed 


as: 





+ YY (=~ + =) ; (23) 


where u and w are the horizontal and vertical components 
of wave-caused veiocity perturbations. The factor K is 
the coefficient of turbulent exchange which is related to 
the generation of turbulent fluctuations by the mean veloc- 
ity gradient. The element V is the coefficient of turbu- 
lent exchange which is related to the wave-caused velocity 
fluctuations. It will be assumed that v >> Me everywhere 
except in the thin laminar sublayer in which we shall ne- 
glect the effect of turbulent viscosity. Both the coeffi- 
cients, K and v, will be assumed to be constant in the 
R-direction and V Jis@a function of height. In jaddition, 
tHe coeh&iicitent K ale be assumed to be given by: 


K = K Uy, Z 


which increases linearly with height and leads to the un- 


disturbed mean velocity profile 


aes a yr Der 
Uae ae eee (2.4) 


r3 





Incorporating the continuity equation with the Bous- 
sinesq approximation, (V*V) = 0, the equations may be 


written in component form as: 


ee eer Ue Boel 8P ieee, 2.) Bs, au 4 fea , aw 

aM Vutwae = Sot + yx ies) -F a7 Kas tv (a + ya) 1» 2-5) 
20 eee D2 ALE ee Ba, ce. 

yo +\WWeVw = 5 + 7x Kap + Vi Goes + 50) I+ yk cae (2.26) 


The vorticity (t) is defined to be: 


2 


w'* but for this work the assumption of: 


Normally a = 2 





will be incorporated into the vorticity equation which may 


now be written in the form: 














a 2 a WV . a @ oma ou — 
ee eer Ona VY tag oe) = 
37 3U 9 3 3° ; 
[Ry tVGy-t ag Vt 2) 
32° 9z2°9z ox? 
Li®the Pcontinuity equation, 
du OW _ 
peeoa 
is used, the vorticity equation becomes: 
zZ 
t ox ae 0x oz 
9° U 3 3 : : 
[kee + v(SE + SS) vet + 28) (2.8) 
9z ox*dz dx? 


The products of the perturbation quantities are neglected 
excCepeswhere they affect the mean flow. This ils consistent 


with the consideration of a single wave number. 
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The mean quantities are balanced separately and sub- 











tracted from (2.8). This results in the two-equation 
system: 
2 2 3 3 
2 4yb py 2B = 2 py Sh + yyy se— 4+ 24), (2,9) 
92° d2z7 Ao Ble ax? 
Oe. , Glebe au 
U 7% + ow a a Gi 77) (2.210) 


Matheele introduction of?stream function notation, the 


velocities can be defined as: 


ec ee ee ov) 


az? ax ° (2.11) 


wrens transtorms (2.9) and (2.10) into; 


F 





dt Vaz Vacx? 7 ES i Qa an Vy Gan 
0 
+ = [v(y_. a eee! = ae ca bs e2 eey) 
We a -_ el ae = hl = Klas 2 13° 


Yefimov (1970) simplified the solution by assuming K and \ 
to be of the same order. This assumption will be incor- 
porated into dqudeton (2.13) by replacing K by v. 

The analytic solution for the special case, U = 0, to 
(2.12) and (2.13) has been considered by Yefimov and 
Pososhkov (1970). Their results have been adapted and in- 
corporated as a test of the numerical technique as outlined 


imeAppendix C. 








The stream function will be defined as: 


b= WoCest)t A(z,t)cos[k(x-ct) ]+ B(z,t)sin[k(x-ct)] (2.14) 


where A(z,t) and B(z,t) are the stream function amplitudes. 
Substitute the expressions for the various derivations 
of w into (2.12) and (2.13) and equate Coefficients Of thic 


cosine terms to produce: 





) 2 3) 
2 -(a,, - Ak?) = (U-c)[Bk> - Bk] + BKUL, + 
| 
Ey oe Wee ee ml (2.15) 
22 ZZ ZZ 


The coefficients of the sine terms are: 





5 2. | 3 
2_(B - BK’) = -(U-c) [Ak™ - A, k] ~ AkUL, + 
2 
aT CoM BE-)) + vile k + tk) CC (2.16) 
2 ZZ ZZ 


The mean velocity equation is: 


k 
Fe 7 oe" 7 Sap i AB oz! . Ue ee eo, 


Equations (2.15) and (2.16) determine the wave-caused 
velocity perturbations and therefore are the basic equa- 
tions for solving the problem. Equation (2.17) is used 


for the calculation of the mean velocity pérturbation u. 
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ig We ae SOLUTION OF THE PROBLEM 


A. INITIALIZATION 
To initiate the solution of equations (2.15) and (2.16), 
the undisturbed mean velocity was assumed to be in the 


form: 


S 


* 


WC) ge alee eS a) (3eay) 


O 
as shown in Figure l. 

The functions A(z,t) and B(z,t) were initialized 
through: 

A(z,0) = A(0,0) sinh [k(H1-z)]/sinh(KH1), C252) 
and B(z,0) = 0 , (353) 
where H is the height of the boundary layer and Hl the 


height one grid point less than H.. 


Bee BOUNDARY CONDITICNS 
Yefimov (1970) defined the boundary conditions for the 


lower boundary of the turbulent layer as: 


u(0,t) akcecos[k(x-ct)], 
w(0,t) = akesin[k(x-ct)]. 


In terms of the stream function the problem was solved 


with the boundary conditions: 


A(O,t) = ac, A, (0, t) = -akc, (3%4) 
B(0O,t) = BL (0,t) = 0, 3,5) 
ACH,t) = A, (H,t) = B(H,t) = BCH, t) =O. C36) 


aby 





Z(t) 





Figure 1. Undisturbed mean velocity profile Ci/ sec). 
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Consequently, the boundary layer was limited to a certain 
height H, where the wave-caused velocity perturbations were 
set equal to zero. 

Yefimov (1970) studied a similar problem in which the 
upper limit of the boundary layer was chosen sufficiently 
high so that it was reasonable to assume the conditions of 
(3.6) to hold. This upper limit was found to be greater 
than one and one-half wavelengths in oerdernsto reach aysta= 


tionary solution. 


C. RESULTANT FACTORS 

As a consequence of the initialization funetions (3.1), 
(3.2) and (3.3); the boundary conditions (3.4), (3.5) and 
(amon, the funetions R@, tye ane B(z.c) Wwerer tine stepped 
forward until a steady state solution was reached. At each 
time step the Reynolds wave stress was computed from: 


“ __:i E 
T= -uw = xk (BLA AB). CSD), 


According to equation (2.17) a non-zero Reynolds stress will 
result in a mean horizontal velocity (u) of the wave-caused 
perturbations. This supplement to the initial undisturbed 


? 


mean velocity may be computed as: 


u=- “ \ tTdz (3.8) 


19 








and applied as a correction to U at each time step by: 


on 


* on 

GS) S ake a 1) + Gee (3.9) 
QO ss 

Dr FINITE DLFFERENCING 


The finite differencing was accomplished through the 


use of centered finite difference schemes for all terms 


except those involving friction. These terms were computed 
at time (t-At). Examples of the schemes are: 
RFS oe. £ ee (3.10) 
Zz 2. Kl K-1‘ ? 
and 
ee ee te LBygil) 
Zo ee K+1 K K-1 


As a result of the finite difference schemes and the 
boundary conditions, the variables at the second grid point 
are specified by: 

GA Ae op oO 
In addition the difference schemes and the boundary condi- 
tions specify that the stream function amplitudes A(z,t) 
and B(z,t) remain constant at both the bottom and top two 


grid points. 


E. NUMERICAL SOLUTION 

Poet ans (2.5) and (2.16) were solved utilizing the 
Prameip le of Rachtmyer (1957) withta forward time step to 
Start the marching procedure. To avoid solution separation, 


the problem was restarted every 50 time steps with a finite 
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difference scheme developed by Matsuno (1966). This was a 
two-step iteration to simulate the Euler backward differ- 
ence method. 

All calculations were performed to double precision 
accuracy and continued until the Reynolds stress (3.7) 


reached a steady state condition which is defined by: 


7 alt 
Ey2 
Zz -5 
Re= >. fz(z,t) — 142,t—1) |} al x Ogee 
Oo 

A space step of 0.25 meters and a time step of 0.0125 
seconds were utilized. Additional accuracy was attained 
in some cases by reducing the space step to 0.125 meters 


with a corresponding time step of 0.00312 seconds. 


Ya 





IV. COMPUTATIONAL PROCEDURES 


A. VARIABLES 

The forecast formulas, as stated in (2.15) and CLP Sp ee 
alow for five independent variables. These are wave 
amplitude (a), wave number (k), dynamic velocity (U,), rough- 
ness length (Zo) and the coefficient of turbulent exchange 
(Vv). 

Yefimov (1970) assumed that the slope of the waves was 
small, so that the expressions for the velocities on the 
boundary predicted by the theory of waves of small slope 
could be applied. This assumption of a small slope is 
applied by establishing the relationship ak = .01, which is 
incorporated in the model as a means of computing the 
initial wave amplitude (3B CON Wie 

The resulting variables k, U,, Zz? and v were then con- 
sidered systematically by first determining the range of 
variation of each of the variables. Then central values 
were selected and used as example cases by holding three 
Wc dabilds fixed and computing test cases of the remaining 
variable over its full range. 

Kraus (1967) studied the relationship of the roughness 
length to that of the mean wind where roughness length is 


given by: 
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The constant, m, was determined by different researchers 
with values ranging from 0.0125 to 0.08. A central value 
of m = 0.042 was selected and the roughness length was com- 


puted from: 


z= 0.00428 Wa ee (4.1) 


for the case where z was dependent upon U,. This reduces 
the independent variables to U,, k and Vv. 
The majority of the experiments were conducted with 


ll 0.0001 meters after Hill (1962). 


B. MODEL CONSIDERATIONS 

Although the value of the coefficient of turbulent ex- 
change is unknown, it is reasonable to assume that it can 
be variable with height as discussed by Yefimov (1970). 
This possibility has been co pora tell in the model and 
several cases have been investigated where V was allowed 


to vary with hetlight. 
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Ven SRESULES 


A. EXPERIMENTAL RESULTS 

The number of variables and their large range of varia- 
tion necessitated numerous numerical experiments. To re-~ 
duce the number of experiments as much as possible, the 
variables were restricted to a limited range. The selection 
of a variation increment was guided by the necessity to 
choose a sufficiently large increment while remaining with- 
im thiessicope of full investigation of the problem. The 
variables and the range of values used are shown in Table l. 
These were used in the various combinations as shown on 
Preuneses = 10. 


TABLE I. Range of the variables considered in the 
numerical experiments. 


*central values 


em) Wy (it /sselc ) on Hise a) Zz fm) 
0.050 O20 0 1.100 Or4s0 
07.07 5 *¥0.020 Sire 4:0 O01 
O00 0.040 0.200 On 00T 
Orai2 5 0.060 0.110 Oo 0005 
O 4155 0 O70S0 0.020 *0.0001 
OR EG AS: *Q0.100 0.011 OF 0000s 
*0,.200 OVA ets OnO00 2 0.00001 
O3225 *0.200 0.000005 
0.2510 0.250 0.000001 
Os 275 0-300 
0.300 0.350 

0.400 





. | 


The existence of a critical level was an important 
factor in the determination of the profile of the solution. 
The critical level is defined as that height where phase 
speed is equal to the mean wind speed. Since the phase 
speed is inversely proportional to the wave number(c = 
AGIA the phase speed will decrease as wave number in- 
creases. Therefore, for a fixed mean wind profile, an 
increase in wave number will decrease the height of the 
critical level (if existing). In the cases of fixed wave 
numbers, an increase in U,, or decrease of Z should de- 
crease the height of the critical level. 

The first case considered was that of a variable wave 
number. In this case the results were greatly affected 
by the value chosen for U, and therefore two experiments 
were conducted. In the first experiment (Figure 2) the 
wave number was varied with U,= 0.20 m/sec. The second 
experiment (Figure 3) was identical but U, = 0.020 m/sec. 

Some noteworthy features of these two tests are: 

a. the Pe ea se in magnitude and variation of T as 
U, was decreased; 

b. the decrease in magnitude and variation of T as 
k was increased; 

Cc. a maximum value of T occurs at approximately one 
meter above the wave surface; 

d. the increase in T at a certain wave number when 


U, = 0.20 m/sec. It is siignificant that this was the 


Zp 








first case in which the critical level was within the layer 
where T varies Significantly. 

The second case considered was that of different mean 
wind profiles. The mean wind given by (3.1) is dependent 
upon -by U, and Zs Since U, and Zz are considered to be 
interrelated there are numerous ways to vary the mean wind 
profile. In this study three experiments were conducted. 
The first case (Figure 4) was that of a fixed wave number, 
roughness length and turbulent coefficient. In the second 
case (Figure 5)the wave number and turbulent coefficient 
remained fixed but the roughness length was allowed to vary 
Pecos to (4.1). Hence, both dynamic velocity and rough- 
ness length varied in this case. The third case (Figure 6) 
was that of fixed dynamic velocity, wave number and turbu- 
lent coefficient, while roughness length was allowed to 
Vary. 

Some notable features are: 

a. the increase in magnitude of the wave-related 
stress as Ze was decreased; 

b. the increase in magnitude of the wave-related 
stress as U, was increased until U, is approximately equal 
BemUa2) m/sec. This is the point where the critical level 
is within the layer where T varies significantly. 

The third case to be considered was that of various 
turbulent coefficients (Figure 7) while the other variables 


were held fixed. 
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A noteworthy feature is the increase in magnitude of T 
as Vv decreased to 0.020 mk eee This is followed by a 
decrease in the magnitude of T as V continued to decrease, 

In the fourth case the turbulent coefficient was 
allowed to vary with height (Figure 8). In the first ex- 
periment the turbulent coefficient was decreased with 
height, to 10 meters, according to the relationship: 

V(z) 9 0.5 tsunh (H-zi)ysinh@). 

On the subsequent experiments V increased linearly with 
height (Figure 8) until v = 0.24 a jade and then remained 
constant thereafter. The results of the first experiment 
are shown in Figure 9 and those of the increasing Vv experi- 
ments are shown in Figure 10. A similar case, but with Vv 
held constant, is shown in Figure 10 for comparison. 

some notable features are: 

a. the similarity in the profiles of the cases of 
increasing V with that of decreasing Vv; 

b. the presence of a maximum at approximately 0.5 
Meters above the wave surface; 

an the rapid change from positive to negative stress; 

d. the decrease in magnitude of T as V increases to 


greater heights. 


B. COMPARISONS WITH DATA 
Additional experiments were conducted to compare the 


results of this model with that of the data analyzed by 


—_ — 





Davidson amd Frank (1973). Th& initial input parametens 
(Table II) were combined with a constant turbulent coeffi- 
cient (v = 0.240 Cy acl) and the calculations time stepped 
forward until a steady state was reached. The resulting 


stress profiles and amplitudes are shown in Figures 11-13. 


TABLE II. Dynamic velocity, wavelength and roughness 
length values determined by Davidson and 
Frank (1973) with corresponding wave numbers 
and critical levels. 


GEteres 


Case HC, sec) \ (m) z (a) Cn) level (m) 
LPs 0.200 35 0.003280 0.179 43.50 
De, 0.167 23 0.000079 0.273 29.00 
3% Oy 124 23 0.000130 0,273 none 


Appendix E outlines the procedure for obtaining data from 
the reported results of Davidson and Frank (1973) and also 
the procedures for interpretation of the numerical resuits. 
The elements for comparison have been displayed in 
mablLe Jit. 

A” significant aeewawrae of Table III is the high-degree 
of correlation of the magnitude of the velocities obtained 
from the variance spectra, phase amplitude and numerical 
amplitude procedures. An equally significant feature is 
the correlation of the magnitude of the stresses obtained 
by the velocity cospectra with those from the numerical 
@alcuratcions. Lt should be noted that the sign of the co- 


spectra stress is dependent upon the wave phase. The 
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differences in the signs should therefore not be considered 
as a Significant discrepancy. 
The comparisons are the first of this kind for data 


from the natural regime. The agreement is significant. 
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Figure 2. Distribution of the wave-caused stress (T in 
m2/sec2) with height for Use 0.2 m/sec, V = 
0.240 m@sec2, z = 0.0001m and k equal, respec- 
tively, to: Seon O> Onmee me Oh O755 3). 100, 
i mO 12 5 oe 0. 150; 6) Om175, 7)8 0.200, 
BO. 225, 9) 0.250, 10)m 0275, 11) 0.2300. The 
arrows indicate the height of the critical level 
and the identifying wave number. 
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Distribution of the wave-caused stress (T in 
m?sec2) with height for U, = 0.02 m/sec, V = 
jemi) sec sO OOUdmmeand ky equal, respec— 
tively, to: oe MMINeEOR ODO. 201.075, 3)0.100, 
Ey Ops? Spano Cre OGG) 0.1755 7). 0.200, 
BOmOI2 2 cS mOnmOr 2 500m 0) 104275, 11) 0.300. 
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Figure 4. Distribution of the wave-caused stress (T in 


n2/sec?) with height for k = 0.200m7l, v = 

0.240 m2/sec, z = '0.0001m, and Us fiquial, nespec- 
tively, to: “SEO. Om ep) 0102 3) 0.04, 

4m O00.) MUR OS. 6) 0 M0, 7) OF le, &) 0.20, 

See tee) pO Oe) LL) OO Joe 2980.40). Bhie 
aELrows! indicate, the helight of the critical level 
and the identifying dynamic velocity. 
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Figure 5. Distribution of the wave-caused stress (T in 
m*/sec2) with height for k = pe? On ie ve—)O. 240 
m2/sec, z = 0.00428 U,* andwue €aual oer espective— 
ven tome mel) 0-01 2pm 0R02, 3) 0704, 4). 0.06; 
SpaOe0 Gum opm 0.10. /)e 0, 15~ 58) 0.20,.9) 0.25, 

LOO sri) O, Jone) On 4G. Thetarrows 
indicate the height of the critical level and 
pices dence may imc dynanie yveLocucy . 
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Figure 6. Distribution of the wave-~caused stress (Tt in 
m2/sec2) with height for U, = 0.100 m/sec, 
k = 0.200 m-1, v = 0.240m2)sec and z_ equal, 
aac Cero meso l) 0e ye2) 0.01, 3!) 0.001, 
4) 0.0005, 5) 0.0001, 6) 0.00005, 7) 0.00001, 
8) 0.000005, 9) 0.000001. 
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Figure 7. Distribution of the wave-caused stress CTaeeien 

m2/sec) with height for k = 0.200 m-4, U, = 
0.100 m/sec, z = 0.0001m, and v equal, 


respectively, “to: 1) 1.100, 2) 0.200, 
ay Omio ea) or020, 5) 0.011, 6) 0.002. 
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Figure 8. 





Variation of the turbulent coefficient with 
height, respectively, to: 1) decreasing, 

2) increasing to 5m, 3) increasing to 10m, 
4) increasing to 15m, 5) increasing to 20m, 
6) increasing to 25m. 
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Pua ee DI Stribueteongea: tCheswave-caused@stress (t in 
m*/sec*) with height for U, = 0.100 m/sec, 
k®=.0.2008n7 4. z = 0.0001m and v decreasing 


with height. 
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Dis tributaomof tthe wave-caused stress (T in 
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height, respectively, Pom Jon om, o2) LOm, 3) 
1500, ) Ole ee 2 Ome Curves Oo mis the case of 
constant v with the same Uae & and Z values. 
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Distribution of the wave-caused stress (T in 
m2/sec2) with height for k = 0.179m71, lie 5 
0.200 m/sec, 2 a 0.003280m, and v = 0.240m*/ 
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Figure 12. Distribution of the wave-caused stress (T in 
MPP oe? ath heueht fork = On273m -, VU, = 
0.167m/sec, Za H2 00007 Sue and Vv = 0.240m2/ 

Seca. 
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Figure 13. 
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Distribution of the wave-caused stress (tT in 
m2/sec2) with height for k = 0.273m 4, U 


- a 
B. pileiee ee B= 0 O00US One anday=0.240 ne 
sec’, e 
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VI. CONCLUSIONS 


The Reynolds wave stress of the perturbation-velocity 
field in the atmospheric near water layer above waves has 
been obtained from a numerical model. The resultant Rey- 
nolds stress can be used to compute the energy which is 
transferred by the average motion of the air to the wave 
fluctuations and vice versa. Thus a positive stress indi- 
cates the wind energy maintains the wave fluctuations sig- 
nifying a transfer of energy from the wind to the waves. 

The maximum Reynolds stress and greatest energy trans- 
fer was found to exist for the largest mean wind profile 
(largest U, and smallest Zs the lowest wave numbers and a 
turbulent coefficient of approximately 0.020 a In 
addition, the maximum stress occurred at approximately one 
meter above the wave surface. 

A striking result is the increase in the magnitude of 
the Reynolds stress as the mean wind velocity approaches 
that of the wave phase speed in the lowest levels of the 
boundary layer. For those cases where the critical level 
occurs in the region of the maximum stress fluctuations 
Emetestsea sipnitieant increase in the magnitude’ of tite 
stress. Thus the energy is greatest when wave velocity is 
close to the mean wind velocity in the lower levels. 

Another significant result is the change which occurred 


in the stress profile when the turbulent coefficient was 
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allowed to vary with height. The net effect of this change 
is a decrease in the total energy transfer. 

The comparison of numerically generated velocities and 
stresses with those of measured data yielded a very high 
degree of correlation. These results will justify future 
numerical research and data collection to provide new in- 
sights into the air-sea energy transfers. Specifically, the - 
correlation of a numerical model with measured wind data 
could determine the characteristics of the turbulent coeffi- 
cient due to the wave-caused motion. Concurrent studies 
being conducted by Aurand (1973) and Stricker (1973) may be 
considered as sources for additional data to correlate the 


Statistically and numerically derived stress patterns. 
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APPENDIX A 


BAROTROPIC INSTABILITY TEST 


Neuman (1969) considered a simple barotropic numerical 
model. The object of the test was to determine the growth 
rate as a function of the wave number. The similarity be- 
tween the present study and that of Neuman allowed for a 
check of the present model for the special case of no 
ted Cd On) 

However, to accomplish the check, several edi ficatmene 
were made to the present model. meee are: 


A(z,0) 


cos es » and 


0. 


vV(z) 


The mean wind velocity profile is given by: 


U(z) = -U, tanh Ci) (Al) 


as shown in figure 14. 

The procedure required marching the calculations for- 
ward in time until the growth rate stabilized. The point. 
of maximum Reynolds stress for three different times were 
then compared. The magnitude of the stress at these points 
was found to be growing exponentially, as was expected. 

As an additional check, the growth rate of the Reynolds 
stress was evaluated and compared to the growth rate as 
given by the function: 


201 


oni (—— =) 2. (A2) 


Yo 
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Figure 14. Mean wind profile for U = -U. tanh(y/y)), 
(m/sec). 
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In this case the stress at the maximum point was equated 


to (A2) and the equation solved for n using: 


U. = 9,0 m/sec , 
ese 4.0m. 
The resultant growth rate was 0.181. This is more 


accurate than that of Neuman who found a growth rate of 
0.160, but it is less than that of Betchov and Criminale 
(1967) who set the rate at n = 0.19. The latter difference 
is because the present solution has been solved with a 


Pinitecec boundary. 
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APPENDIX B 


EXAMINATION OF THE VISCOSITY TERMS 


In this test the initial field of the cosine stream 
function amplitudes was computed utilizing equation (3.2). 
At this point the mean wind and phase speed were set to 
zero. Then the forecast equations were marched forward 
until a steady state was reached. 

These numerical values were then compared with the 
analytical values derived from the solution of the equa- 


tion: 


: 2A + Ag 0 . (B1) 
ZAZA, 72, ZZ 


The boundary conditions utilized in the solution of (B1) 
are similar to those given by (3.4) and (3.6). These 


Conditions $ane: 


A(0) 


Wig 6 A, (0) = -2ke , 


I 


A(H) A, (CH) = 0. 


The solution of (Bl) is given by: 


hz ir — i ceniscimckz jet cc. ecosi( kz) + (¢ 


Z 3 


1 kz ean le) 


+ kz cos(kz) , (Bz) 


ih 
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where: 


c, = = [A (0) - ¢, kl , 
ys h 
C4 = - FE Sin THE) (A, (0) sin (Hk) + A(0)cos(Hk) + 
c, sin(Hk)+ c, Hk cos (Hk)] , 
c, = i (icon ve hein (atk) 
+ (AS + a(0)k cot (Hk) I cos (Hk) }] 


The comparison of the numerical results with the 
analytical results derived from (B2) indicated that the 
Pininon occurred at the same height for both cases. How- 
ever, the magnitude of the numerical minimum was smaller 


than that of the analytical results. 
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PAPE RIND IT Xs C 


EXAMINATION OF THE CASE OF NO MEAN WIND 


Yefimov and Pososhkov (1970) derived an analytic solu- 
tion to the forecast equations (2.15) and (2.16) for the 


Special case U = 0. This solution is given by: 





tT = -u'w' = a ae + ae 2K/ £0 cos (SS : (C1) 
20 
where: 
Vk 
OS Ss e 
c 
The condition U = 0 was established and the numerical 


solutions of (2.15) and (2.16) were computed to a steady 
state condition. The results of the Priaily bacail and numeri- 
cal calculations of stress were compared for several dif- 
ferent Wave numbers and for different values of the turbu- 
lent exchange coefficient. 

The results of the stress obtained by both methods 
were very Similar. In all cases compared the analytical 
stress reached a minimum at a lower height and somewhat 
smaller magnitude than the numerically derived stress. 
This discrepancy, however, has been primarily attributed 


to the scale of the grid increment. 
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APPENDIX D 


EXAMINATION OF THE RICHTMYER SCHEME 


This test was devised solely to test the method of 
solving equations as given by Richtmyer (1957). The test 


was conducted by equating the right hand side terms of 





equation (2.15) to Sin (=) where w = (H-3)Az. The result 
is: 
Z 
(2 - &?) 24 = sin (FY , 
92:7 S he 


and upon solution yields: 


a ee ae Die | (D1) 


The numerical solution was marched forward one complete 
time step. Then the minimum value of the stream function 
amplitude, A(z,t), was compared with the value derived 
from (D1). The result was a comparison of the values 


LOewe thine thee @pier cemt@ote thle Moiiid erences, 
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APPENDIX E 


EVALUATION OF ANALYZED DATA 


Davidson and Frank (1973) reported observed values of 
Wave heights and wel oct ee fluctuations at two levels above 
the waves. The results from their spectral and joint pro- 
bability density function, conditional mean function (JPDF- 
CMF) analyses are shown in Figures 15 and 16. 

The velocity variance spectra graphs show plots of 
n *° @?(n) vs. In(n) where n is the frequency of the wave in 
Hz. The area under the significant spikes of the curve 
(shaded areas of Figure 15) can be used to determine the 
value of the square of the wave-caused velocity fluctua- 


[6 ab (elal qa or a): These areas may be computed by: 


Ge mettre OIC Tl ) ery TLC TL (El) 


where A In(n) is defined at the significant peak. 

The velocity cospectra graphs (Figure 15) are similar 
to those above but the area in this case is indicative of 
the wave-caused Reynolds stress (u'w'). These results 
have been presented in Table III for comparison. 

Davidson and Frank (1973) also presented phase-~ampli- 
tude results for the same data as above. These results 
(Figure 16), which may be interpreted directly, are shown 


in, Tab Mey LEL for compariison. 


52 





The numerical maximum amplitude has been computed 
from the steady state values of the stream function ampli- 


tude by: 


' 2 B 21/2 
z Zz 


= 
il 
a 
> 
= 


(E2) 


and 


w' k (A? + By (E3) 
which follow as a result from equations (2.11) and (2.14). 
The numerical results were computed with the small wave 
slope assumption where ak = 0.1. However, this exact con- 
dition is not normally met in nature. Therefore, the 
values obtained from equations (E2) and (E3), as well as 
the numerically derived Reynolds stress, must be scaled 
to reflect the actual wave slope. 

The calculation of a scaling factor first involved the 
determination of the actual amplitude (ay)- This may be 
interpreted directly from the "waves" diagram in Figure 16. 


The amplitude used for the numerical calculations may be 


computed by: 


The resulting velocity scaling factor (Fs) may now be ex- 
pressed by: 


a 
Fs = 2 ; (E4) 
N 


Fq = -—> ; Cao) 
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